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Abstract 

We analyze stability of a fermion system with model repulsive pair 
interaction potential. The possibility for different types of restructur- 
ing of the Fermi ground state (at sufficiently great coupling constant) 
is related to the analytic properties of such potential. In particular, for 
the screened Coulomb law it is shown that the restructuring cannot be 
of the Fermi condensation type, known earlier for some exactly solvable 
models, and instead it belongs to the class of topological transitions 
(TT). For this model, a phase diagram has been built in the variables 
"screening parameter - coupling constant" which displays two kinds of 
TT: a 5 /2-kind similar to the known Lifshitz transitions in metals, and 
a 2-kind characteristic for a uniform strongly interacting system. 

PACS: 71.10.-w, 71.10.Hf 

The common ground state of isotropic Fermi gas with density p is de- 
scribed (at zero temperature) by the stepwise Fermi function n-p(p) = 9(p-p — 
p), dropping discontinuously from 1 to at the Fermi momentum p-p = 
(3ti 2 p) 1 ^ 3 (in the units where % = 1). The Landau theory of an interacting 
Fermi liquid started from the assumptions that the quasiparticle distribution 
function n(p) coincides with the one of an ideal gas, while the single particle 
spectrum e p is similar to that of ideal gas, being characterized by the effective 
mass M* gj. 



But these and e p can be broken down under certain circumstances. 

The best known example is the Cooper pairing at arbitrarily weak attractive 
interaction with subsequent formation of the pair condensate and gapped 
quasiparticle spectrum [|J, but also a repulsive interaction, if sufficiently 
strong, can produce non-trivial ground states. The first example of such re- 
structuring for a Fermi system with model repulsive interaction [|] revealed 
existence of a critical value a cr for the interaction constant a, such that at 
a = a cr the stability criterion s(p) = (e p — e-p)/(p 2 — Pp) > fails just at 
the Fermi surface: s(pp) = (pp-instability) . Then at a > a cr an exact 
solution exists of a variational equation for n(p) (following from the Landau 
energy functional E[n(p)]), exhibiting a certain finite interval (pi,P2) around 
Pp where the distribution function n(p) varies continuously and takes inter- 
mediate values between 1 and 0, while the single-particle excitation spectrum 
e p has a flat plateau, 
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with (jl being the chemical potential. It is seen from eq. ([I]) that the occupa- 
tion numbers n(p) become variational parameters, deviating from the Fermi 
function in order to minimize the energy E. This phenomenon was called 
fermion condensation (FC) and soon several model forms for E[n(p)} were 
proposed providing similar solutions. However, these models do not 

answer the question whether there exist other types of phase transitions re- 
lated to a rearrangement of the Fermi function np(p) and the single particle 
spectrum e p . It is pertinent to note that the idea of multiconnected Fermi 
sphere, with production of new, interior segments of the Fermi surface, has 
been considered already || pT0[| . This determines an interest to examine the 



overcritical regimes a > a cr for the models displaying alternative types of 
instability. 

The main goal of this Letter is to consider possible types of rearrangement 
of the Fermi ground state. We shall relate the types of such rearrangement 
to analytic properties of the single particle potential of a system in question 
and show that there can exist phase transitions of different type from FC. 

The general scheme of Refs. for a homogeneous system of non- 

polarized fermions with mass M and model isotropic interaction potential 
U(p) considers the energy functional E[n( 
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EHP)} = /^«W(^J + i//nW^(IP "Pi)-*') f^. (2) 
and the related quasiparticle dispersion law: 

Performing the angular integration and passing to the dimensionless vari- 
ables: x = p/p-p, V — ZMep/pp, z = 2ti 2 ME/py, eqs. (||,|) can be presented 
in a simpler form: 

z[v(xj\ = y [x 4 + ]-x 2 V(x)\v(x)dx, (4) 

y(x) = X 2 + V(x), (5) 

with the single particle potential V(x) being given by, 

V(x) = —J x'v(x')u(x,x')dx', (6) 

, , x+x' 

u(x,x) = — — / ultUdt. 
n 2 p F J 

\x-x'\ 

Here u(x) = [/(ppx), and the distribution function u(x) = n{p-px) is positive, 
obeys the normalization condition: 

x 2 v{x)dx = 1/3, (7) 

and the Pauli principle limitation v(x) < 1. The latter can be lifted using, 
e.g., the ansatz: v(x) = [1 + tanhry(x)]/2, then the system ground state 
should correspond to the minimum of the functional: 



f[v( x )] — J [1 + tanh?7(x)]{x 4 — [ix 2 + -x J x'[l +t&nhr}(x')]u(x,x')dx'}dx, 

(8) 
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containing a Lagrange multiplier /i, with respect to an arbitrary variation of 
the auxiliary function r)(x). This permits to present the necessary condition 
of extremum 5f = as: 

x 2 u(x)[l-u(x)][y(x)-fi] = 0, (9) 

which means that either u(x) takes only the values and 1 or the dispersion 
law is flat 0: y(x) = /i, in accordance with eq. (|l|). The last possibility just 
corresponds to FC. As it is seen from eq. (|l|), the single particle spectrum 
e p in this case cannot be analytic function of complex p in any open domain, 
containing the FC interval Ap = [pi,P2]- In fact, all the derivatives of e p 
with respect to p should be zero along Ap, while this is not the case along 
the real axis outside Ap. For instance, in the FC model with U(p) = U /p 
the kernel, eq. (|]), results non-analytic: 

u(x, x ) = — — (x + X + \x — X \) (10) 

which eventually causes non-analyticity of the potential V(x). 

On the other hand, as it follows from eq. (H), the single particle spectrum 
will be an analytic function along the whole real axis, provided V(x) is such 
a function. In this case FC is forbidden and the only alternative to the Fermi 
ground state (if the stability criterion gets broken) leaves in a topological 
transition (TT) between the topologically unequal states with v(x) = 0, 1 
JT2|. Generally, all such states are classified by the indices of connectedness 
(known as Betti numbers in algebraic topology |13|) for the support of u(x). 
In fact, for an isotropic system, these numbers are simply to count the sep- 
arate (concentric) segments of the Fermi surface. Then the system ground 
state will correspond to such a multiconnected distribution (Fig. 1): 

n 

u ( x ) = ~ x 2l -i)9(x 2i - x), (11) 

i=i 

that the parameters < X\ < x 2 < ■ ■ ■ < x 2n obey the normalization 
condition: 
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Figure 1: Occupation function for a multiconnected distribution. 



and the related z, eq. (|4|), 

z= 2^ J A^ + V^dx (13) 

has the absolute minimum with respect to x\, . . . , x^n-i and to n > 1. 

Differentiating eq. ( |I3| ) in the parameters x%, . . . , X2 n -i with use of the 
relations: 

= l<Jfe<2n-l, (14) 

by virtue of eq. (|T2|), and taking into account that the potential V(x) in the 
dispersion law y(x) in fact also depends on these parameters: 



1 ^ /'•2'2i 

V( x ) = -J2 x'u(x,x')dx', (15) 

3^ J X2i-1 



we present the necessary conditions of extremum as: 



-l) fc x^[y(x fe ) - y(x 2n )] = 0, 1 < A; < 2n - 1. (16) 



Hence the multiconnected ground state is controlled by the evident rule of 
unique Fermi level: y(xk) = y(x2 n ), for all 1 < k < 2n — 1 (except for x± = 0). 
In principle, given the kernel u(x, x') all the (2n — 1) unknown parameters Xk 
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can be found from eqs. (|16|,|5|,|T5|) . Then, the sufficient stability conditions: 
d 2 z / dxidxj = ctidij, «j > 0, provide the generalized stability criterion: the 
dimensionless function 

„ W = 2M.fr) = (17) 

should be positive within filled and negative within empty intervals, turning 
zero at their boundaries in accordance with eqs. fll6|). It can be proved 
rigorously that for given analytic u(x, x') this criterion uniquely defines the 
system ground state. 

In what follows we shall label each multiconnected state, eq. (|TT|), by an 
entire number related to the binary sequence of empty and filled intervals 
read from x 2n to 0. Thus, the Fermi state with a single filled interval [x 2 = 
1, X\ = 0] reads as unity, the state with a void at origin: (filled [x 2 , Xi] and 
empty [x\, 0]) reads as (10) = 2, the state with a single gap: (101) = 3, etc. 
Note that all even phases have a void at the origin and odd phases have not. 

For free fermions V(x) = 0, y(x) = x 2 , eqs. (|16D only yield the trivial 
solution corresponding to the Fermi state 1. In order to pass to non-trivial 
realizations of TT, we choose U(p) corresponding to the common screened 
Coulomb potential: 

47TP 2 

u(p) = -r— 2 . (is) 

The related explicit form for the kernel 

< X ^ =a ^ %V^)2 + + % ( 19 ) 

with the dimensionless screening parameter Xq = Po/p-p and the coupling 
constant a = 2Me 2 /iip-p, evidently displays the necessary analytical proper- 
ties for existence of TT. 

With use of eqs. (p"5|) and (0), the potential V (x) is expressed in ele- 
mentary functions: 



i=i 



V(x;x',x") = a{2 



V{x)=Y,V{x]X2i-l,X2i), (20) 

2x (x" — x') (x 2 + Xq + x'x") 



x" — x' — xq arctan 



(X 2 + Xq — X' 2 ) (X 2 + Xq — X" 2 ) 



+ 
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+x 2_ x 2\ ln (X + X'f+Xl _ f, 2 - ^ (X + X'f+Xl 

(X — x") + Xn [X — X') + Xqj 



1 

+ 2x~ 

Then, the straightforward analysis of eqs. (|16[) shows that their non-trivial 
solutions only appear when the coupling parameter a exceeds a certain criti- 
cal value a*. This corresponds to the moment when the stability criterion || 
a(x) = (yp(l) — y-p(x))/(l — x 2 ) > calculated with the Fermi distribution, 
y-p(x) = x 2 + V(x; 0, 1), fails in a certain point < Xi < 1 within the Fermi 
sphere: a(xi) — > 0. There are two different types of such instability depend- 
ing on the screening parameter x (Fig. 2). For xo below certain threshold 
value x t h ~ 0.32365 (weak screening regime, WSR) the instability point Xi 
sets rather close to the Fermi surface: 1 — Xj -C 1, while it drops in a critical 
way to zero at xq — > x t h and pertains zero for all xq > x t h (strong screening 
regime, SSR). The critical coupling a*(xo) results a monotonically growing 
function of Xo, having the asymptotics o*{xq) ~ (ln2/a;o — l) _1 at xq — > 
and staying analytic at a t h = o>*(x t h) ~ 0.91535, where it only exhibits an 
inflexion point. 

These two types of instability give rise to different types of TT from the 
state 1 at a > a*: at SSR a void appears around x = (1 — > 2 transition), 
and at WSR a gap opens around Xi (1 — ► 3 transition). Further analysis 
of eqs. flTB]) shows that the point x t h, octh represents a triple point in the 
phase diagram in variables x , a (Fig. 3) where the phases 1, 2, and 3 
meet. Similarly to the onset of instability in the Fermi state 1, each TT to 
higher order phases with growing a is manifested by that cr(x), eq. (|17D, 
turns zero at some point < Xi < Xm different from the existing interfaces. 
If this occurs at the very origin, X{ = 0, the phase number rises at TT 
by 1, corresponding to opening of a void (passing from odd to even phase) 
or to emerging "island" (even — > odd). For Xi > 0, a thin spherical gap 
opens within a filled region or a thin filled spherical sheet emerges within 
a gap, then the phase number rises by 2, not changing the parity. A part 
of the whole diagram shown in Fig. 3 demonstrates that with decreasing 
Xo (weaken screening) all the even phases terminate at certain triple points. 
This agrees in particular with the numeric study of the considered model 
along the line x = 0.07 at growing a ||10| , where only the sequence of odd 
phases 1 — ► 3 — > 5 — > . . . was indicated (shown by the dashed arrow in Fig. 
3). 
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Figure 2: Critical coupling a* and instability point Xi as functions of screen- 
ing. The regions of weak screening (WSR) and strong screening (SSR) are 
separated by the threshold value x t h (note that x t h, &th is the triple point 
between the phases 1, 2, 3 in Fig. 3). 
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Figure 3: Phase diagram in variables "screening-coupling". Each phase with 
certain topology is labeled by the total number of filled and empty regions 
(see Fig. |l|). Even phases (shadowed) are separated from odd ones by " 5 /2- 
kind" topological transition (TT) lines, while 2-kind TT lines separate odd 
phases. Triple points, where two 5 /2-TT and one 2-TT meet, are shown by 
circles. The dashed arrow indicates the phase trajectory studied numerically 
in Ref. 0. 
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The energy gain A (r) at TT as a function of small overcriticality param- 
eter r = a/a* — 1 is evidently proportional to r times the volume of a new 
emerging phase region (empty or filled). The analysis shows that the radius 
of a void (or an island) for a TT with change of parity is ~ t 1 / 2 . Conse- 
quently we have A (r) ~ r 5//2 which indicates a similarity of this situation 
to the known " 5 / 2 -kind" phase transitions in the theory of metals |fPH , but 
its specific character is that the new segment of Fermi surface opens at very 
small momentum values, which can dramatically change the system response 
to, e.g., electron-phonon interaction. On the other hand, this segment may 
have a pronounced effect on the thermodynamical properties of 3 He at low 
temperatures, especially in the case of P-pairing. 

For a TT with unchanged parity, the width of a gap (or a sheet) is found 
~ r, hence the energy gain results A (t) ~ t 2 , and such TT can be assigned 
2 kind. It follows from the above consideration that each triple point in the 
x — a phase diagram is a point of confluence of two 5 / 2 -kind TT lines into one 



2-kind line. The latter type of TT was already discussed in literature [y, |10 
and we only mention here that its ocurrance on a whole continuous surface 
in the momentum space is rather specific for systems with strong fermion- 
fermion interaction, while the known TT's in metals, under the effects of 
crystalline field, occur typically at separate points in the quasimomentum 
space. 

It is of interest to note that in the limit xq — > 0, a — > 0, reached along 
a line a = kxo, we attain the exactly solvable model: U(p) — > (2ii) 3 Uo5(p) 
with U = fc/(2Mp F ), which is known to display FC at all U > |J. The 
analytic mechanism of this behavior consists in that the poles of U(p), eq. 
(|18| ) tend to zero, thus restoring the analytical properties necessary for FC. 
Otherwise, the FC regime corresponds to the phase order — > oo, when the 
density of infinitely thin filled (separated by empty) regions approaches some 
continuous function < v(x) < 1 and the dispersion law turns flat by 
eq. (|1|). 

A few remarks are in order at this point. First, the considered model 
formally treats xq and a as independent parameters, though in fact a certain 
relation between them can be imposed. Under such restriction, the system 
ground state should depend on a single parameter, say the particle density p, 
along a certain trajectory a(x ) in the above suggested phase diagram. For 
instance, with the simplest Thomas-Fermi relation for the free electron gas: 
a(xo) = xl/2, this trajectory stays fully within the Fermi state 1 over all the 
physically reasonable range of densities. Hence a faster growth of a(xo) is 
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necessary for realization of TT in any fermionic system with the interaction, 
eq. flU). 

Second, the single particle potential V(p) of a real system cannot be an 
entire function of p around p-p because of the stepwise form of the quasipar- 
ticle distribution. Therefore, as the coupling constant moves away from the 
critical value a* within the WSR domain, the concentric Fermi spheres will 
be taken up by FC. A close look at the role of the density wave instabil- 
ity, which sets in at sufficiently large a, shows that this is true ||. In fact, 
these arguments do not work in the case of SSR. Thus, it is quite possible 
to observe the two separate Fermi sphere regimes. There is a good reason to 
mention that neither in the case when the FC phase transition takes place nor 
in the case when types of TT are present the standard Kohn-Sham scheme 



|Tj| is no longer valid. Beyond the FC or TT phase transitions the occu- 
pations numbers of quasiparticles serve as variational parameters. Thus, to 
get a reasonable description of the system, one has to consider the ground 
state energy as a functional of the occupation numbers E[(n(p)] rather then 
a functional of the density E[p] ||T5| . A more detailed study of such systems, 
including the finite temperature effects, is in order. 
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